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Abstract. We compute the luminosity function (LF) and the formation rate of long gamma ray 
bursts (GRBs) in three different scenarios: i) GRBs follow the cosmic star formation and their 
LF is constant in time; ii) GRBs follow the cosmic star formation but the LF varies with redshift; 
iii) GRBs form preferentially in low-metallicity environments. We then test model predictions 
against the Swift 3-year data, showing that scenario i) is robustly ruled out. Moreover, we show 
that the number of bright GRBs detected by Swift suggests that GRBs should have experienced 
some sort of luminosity evolution with redshift, being more luminous in the past. Finally we 
propose to use the observations of the afterglow spectrum of GRBs at z ^ 5.5 to constrain the 
reionization history and we applied our method to the case of GRB 050904. 
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1. Introduction 

Long Gamma Ray Bursts (GRBs) are powerful flashes of high-energy photons occur- 
ring at an average rate of a few per day throughout the Universe up to very high redshift 
(the current record is z = 6.29). The energy source of a long GRB is believed to be 
associated to the collapse of the core of a massive star (see Meszaros 2006 for a review). 
One of the main goals of the Swift satellite (Gehrels et al. 2004) is to tackle the key issue 
of the GRB luminosity function (LF). Unfortunately, although the number of GRBs with 
good redshift determination has been largely increased by Swift, the sample is still too 
poor (and bias dominated) to allow a direct measurement of the LF. We use here the 
Swift 3-year data to constrain the GRB LF and its evolution (Salvaterra & Chincarini 
2007, Salvaterra et al. 2008b). Moreover, we show a possible use of GRBs detected at 
z ^ 5.5 to study the history of reionization (Gallerani et al. 2008). 



2. Model description 

The observed photon flux, P, in the energy band E m i n < E < E max , emitted by an 
isotropically radiating source at redshift z is 

n _ ^+^sS:^::s(E)dE 

where S{E) is the differential rest-frame photon luminosity of the source, and c?l(z) is the 

119 



120 



R. Salvaterra et al. 



luminosity distance. To describe the typical burst spectrum we adopt the functional form 
proposed by Band et al. (1993), i.e. a broken power-law with a low-energy spectral index 
a, a high-energy spectral index /3, and a break energy — (a — [3)E p /(2 + a), with a = 
— 1 and (3 = —2.25 (Preece et al. 2000). In order to broadly estimate the peak energy of the 
spectrum, E p , for a given isotropic-equivalent peak luminosity, L = J^^^^ ES(E)dE, 
we assumed the validity of the correlation between E p and L ( Yonetoku et al. 2004) . 
Given a normalized GRB LF, (j)(L), the observed rate of bursts with Pi < P < Pi is 

at J dz 4?r 1 + z Jl(p uZ ) 

where dV(z)/dz is the comoving volume element[j], Af2 s is the solid angle covered on 
the sky by the survey, and the factor (1 + z)^ 1 accounts for cosmological time dilation. 
Finally, *grb(z) is the comoving burst formation rate. In this work, we assume that the 
GRB LF is described by a power law with an exponential cut-off at low luminosities, i.e. 
cj>(L) oc (L/L cut )~Z exp(-L cut /L). 

We consider three different scenarios: i) no evolution model, where GRBs follow 
the cosmic star formation and their LF is constant in time; ii) luminosity evolution 
model, where GRBs follow the cosmic star formation but the LF varies with redshift; 
iii) density evolution model, where GRBs form preferentially in low-metallicity en- 
vironments. In the first two cases, the GRB formation rate is simply proportional to the 
global SFR, i.e. ^grb(z) — /cgrb*^(^)- We use here the recent determination of the SFR 
obtained by Hopkins & Beacom (2006), slightly modified to match the observed decline 
of the SFR with (1 + z)~ 3 - 3 at z > 5 suggested by recent deep-field data (Stark et al. 
2006). For the luminosity evolution model, we also assume that the cut-off luminosity in 
the GRB LF varies as L cut = L {1 + z) s . Finally, for density evolution case, the GRB 
formation rate is obtained by convolving the observed SFR with the fraction of galaxies 
at redshift z with metallicity below Z t h using the expression computed by Langer & 
Norman (2006). In this scenario, L cut — const = Lq. 



3. Results 

The free parameters in our model are the GRB formation efficiency fcoRB, the cut-off 
luminosity at z — 0, Lo, and the LF power index £. We optimized the value of these 
parameters by x 2 minimization over the observed differential number counts in the 50- 
300 keV band of BATSE (Stern et al. 2000). We find that it is always possible to find a 
good agreement between models and data. Moreover, we can reproduce also the 3-year 
differential peak flux count distribution in the 15-150 keV Swift band without changing 
the best fit parameters (Salvaterra & Chincarini 2007). We then check the resulting 
redshift distributions in the light of the Swift 3-year data, focusing on the large sample 
of GRBs detected at z > 2.5 and z > 3.5 (Fig. 1 panels a & b). The no evolution model is 
ruled out by the number of sure high-z GRBs. This result is robust since does not depend 
on the assumed SFR at high-z nor on the faint-end of the GRB LF. In conclusion, the 
existence of a large sample of bursts at z > 2.5 in the Swift 3-year data imply that GRBs 
have experienced some kind of evolution, being more luminous or more common in the 
past (Salvaterra & Chincarini 2007). 

In order to discriminate between luminosity and density evolution models, we compute 
the number of luminous GRBs, i.e. bursts with isotropic peak luminosity L ^ 10 53 erg 

f We adopted the 'concordance' model values for the cosmological parameters: h — 0.7, 
Q, m = 0.3, and D. A = 0.7. 
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s -1 in the 1-10000 keV band (Salvaterra et al. 2008b). We compare model predictions 
with the number of bright bursts detected by Swift. Conservatively, our data sample 
contains only bursts with a good redshift measurement and whose peak energy was well 
constrained by Swift itself or other satellites (such as HETE-2 or Konus-Wind). We 
stress here that this number represents a lower limit on the real number of bright GRBs 
detected, since some luminous bursts without z and/or E p can be present in the Swift 
catalog. Results for the pure luminosity (density) evolution models are plotted in the 
panel c (d) of Fig. 1. Data are shown with the histogram where the shaded area takes into 
account errors on the determination of L. We find that models involving pure luminosity 
evolution requires <5 > 1.5 to reproduce the number of known bright GRBs. On the other 
hand, models in which GRB formation is confined in low-metallicity environments fall 
short to account for the observed bright GRBs for Zth > 0.1. Assuming Zth = 0.1 Z©, 
the model reproduces the observed number of bright GRBs, taking also into account the 
errors in the determination of L. This means that essentially all bright bursts present 
in the 3-year Swift catalog have a measured redshift and well constrained peak energy. 
So, although this model can not be discarded with high confidence, the available data 
indicate the need of some evolution in the GRB LF even for such a low value of Z t h- 
For Z t h = 0.3 Z Q , as required by collapsar models (MacFadyen & Woosley 1999), only 
~ 6 bursts with L 10 53 erg s _1 should have been detected in three year, largely 
underpredicting the number of Swift sure identifications. Thus, pure density evolution 
models, where the GRB LF is constant with redshift, are ruled out by the number of 




Figure 1. Panels a & b: cumulative number of GRBs at z > 2.5 (a) and at z > 3.5 (b) 
as a function of the photon flux P. Dotted line refers to the no evolution model, short dashed 
to the luminosity evolution model (S = 1.5) and long-dashed to the density evolution model 
(Zth = 0.1 Zq). The number of sources detected by Swift in three years is shown as solid 
histogram. Note that the observed detections are lower limits, since many high-z GRBs can be 
missed by optical follow-up searches. A field of view of 1.4 sr for Swift is adopted. Panels c 
8z d: cumulative number of luminous GRBs detected by Swift in three years, shown with the 
histogram, as function of the isotropic equivalent peak luminosity, L. Shaded area takes into 
account the errors on the determination of L. Note that the data are to be considered as lower 
limits of the real number of Swift detections. For pure luminosity evolution models (panel c): 
solid line is for 5 = 3, dotted line for 5 = 2.5, short-dashed line for 5 = 2, and long-dashed for 
5 — 1.5. For pure density evolution models (panel d): Solid line is for Zth = 0.1 Zq, dotted line 
is for Zth = 0.2 Zq, and short-dashed line is for Z t h = 0.3 Zq. 
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bright GRBs. In conclusion, available data suggest that GRBs have experienced some 
luminosity evolution with cosmic time. 

4. GRBs from the reionization epoch 

We can now compute a robust lower limit on the number of bursts detectable by Swift 
at very high-z. Assuming a trigger threshold P ^ 0.4 ph s _1 cm~ 2 , at least ~ 5 — 10% 
of detected GRBs should lie at z ^ 5, with > 1 — 3 GRB yr _1 at z ^ 6. These numbers 
double by lowering the Swift trigger threshold by a factor of two (Salvaterra et al. 2008a). 

High-z GRBs are a useful and unique tool to study the Universe near and beyond 
the reionization epoch. Gallerani et al. (2008) have studied the possibility to constrain 
the reionization history using the statistics of the dark portions (gaps) produced by 
intervening neutral hydrogen along the line of sight (LOS) in the afterglow spectra of 
GRB at z ^ 5.5. Two reionization models, both consistent with available observations of 
the high-z Universe, are considered: (i) early reionization model (ERM) where z re ion ~ 7 
and (ii) late reionization model (LRM) where z re i on ~ 6. Suppose now that a GRB at 
rcdshift zqrb is observed at a given flux level in the J band, Fj. We can then ask what is 
the probability that the largest of the dark gaps in its afterglow spectrum is found within 
a given width range. The results are shown in Fig. 2 for two different width ranges; the left 
(right) panels refer to the ERM (LRM) case. The isocontours correspond to a probability 
of 15%, 30%, 45%, and 60%. We find that the two models populate the (zgrb, -Fj) plane 
in a very different way. In particular, for largest gaps in the 40-80 A range, the highest 
probability is obtained for fainter afterglows in the ERM than for the LRM. For largest 
gaps in the range 80-120 A, the probability is in general higher in the LRM with respect 
to the ERM. Note that, in the ERM, only a few spectra should contain the largest gap in 
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Figure 2. Isocontours of the probability that the afterglow spectrum of J-band flux Fj associated 
with a GRB at redshift zgrb, contains the largest gap in the range 40-80 A (top panels) and in 
the range 80-120 A (bottom panel). The left (right) panel shows the results for the ERM (LRM). 
The isocontours correspond to probability of 15% (long dashed line), 30% (short dashed line), 
45% (dotted line), and 60% (solid line). The black point indicates the position in the (zgrb, Fj) 
plane of GRB 050904. 
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this range for Fj > 10 — 40 ^Jy. Fig. 2 allows a straightforward comparison between data 
and model results. It is then natural to apply this procedure to GRB 050904 (black filled 
circle in Fig. 2). The probability to find the largest gap of 65 A is > 45% in the ERM, i.e. 
almost half of the LOS contains the largest gap in the range 40-80 A for a burst with the 
redshift and flux of GRB 050904. Such probability drops for the LRM to ~ 15% clearly 
indicating that in this case the GRB 050904 observation represents a much rarer event. 
Although a large sample of high-z GRBs is required before we conclude that a model 
in which reionization was complete at z ~ 7 is favored by the data, the discriminating 
power of the proposed method is already apparent. 

This kind of analysis requires high signal-to-noise, high resolution spectra of GRB 
afterglow spectra at z ^ 5.5 obtained with the largest ground telescopes soon after the 
burst detection. To avoid wasting observing time, we developed a very effective strategy 
to spot reliable z ^ 5 candidates on the basis of promptly available information provided 
by Swift (Campana et al. 2007, Salvaterra et al. 2007). The selection criteria adopted 
are: long burst observed durations (Tgo > 60 s), faint 7-ray photon fluxes (P < 1 ph s _1 
cm -2 ), and no optical counterpart in the V and bluer filters of UVOT (V > 20). We 
tested our selection procedure against the last ~ 2 years of Swift data showing that our 
method is very efficient and clean (i.e. no low- z interloper is present in the sample). 

5. Conclusions 

We have tested different formation and evolution scenarios for long GRB against the 
3-year Swift dataset. We found that Swift data strongly rule out models in which GRBs 
follow the cosmic star formation and their LF is constant in time. In particular, the num- 
ber of bright GRBs suggests that GRBs should have experienced some sort of luminosity 
evolution with cosmic time, being more luminous in the past. Finally we have shown that 
GRBs at z ^ 5.5 can be use to constrain the reionization history and we applied our 
method to the case of GRB 050904 at z = 6.29. 
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